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Abstract. In this paper we study the filtration laws for the polymeric flow in a porous medium. We use 
the quasi-Newtonian models with share dependent viscosity obeying the power-law and the Carreau's law. Using 
the method of homogenization in [^| the coupled micro-macro homogenized law, governing the quasi-newtonian 
flow in a periodic model of a porous medium, was found. We decouple that law separating the micro from the 
macro scale. We write the macroscopic filtration law in the form of non-linear Darcy's law and we prove that the 
obtained law is well posed. We give the analytical as well as the numerical study of our model. 



1 Introduction 

It is well-known that the viscosity of a polymer melt or a polymer solution significantly changes with the 
share rate (see e.g. jjj). Therefore the Newtonian model, characterized by a constant viscosity, is not a 
reasonable choice for modelling the polymeric flow. The simplest way to overcome that difficulty is to 
use the notion of quasi-Newtonian flow, i.e. to modify the Newton's model by using the variable viscosity 
that depends on the share-rate according to some empirical law. In this paper, in order to describe the 
polymer flow, we use the share dependent viscosities obeying the power-law 

r)(e(u)) = p\e{u)\ r ~ 2 e{u) , 1< r < 2 , p > (1) 

or the Carreau's law 

T)(e(u)) = (r; - r /oo )(l + A|e(w)| 2 ) r/2 - 1 + Voo , 1< r < 2 , 770 > Voo > , A > , (2) 

where 

e(u) = — (Vu + Vu*) the rate-of-strain tensor , 

n 

\ e ( u ) \ = ( \ e ( u )ij | 2 ) 1/ ' 2 — the share rate , 

To find the effective law describing the polymer flow through a porous medium f2 e [] Bourgeat and Mikelic 
in H (see also || and E^] ) used the porous medium with periodic geometry and the homogenization 
technique called the two-scale convergence. Starting from the microscopic problem 

-2div{?7(e(u e ))e(u e )} + Vp £ = e 7 / in fl E (3) 
divu £ = in tt e (4) 
if = on dVl e , (5) 

depending on 7, the following three types of homogenized problems were rigorously derived: 

1. If the Carreau's law was taken for the viscosity and 7^1, the homogenized law is the Darcy's law 

v(x)=K(f(x)-\7 x p°(x)) , ie!! (6) 
div x v = in Q , (7) 

where is the whole domain (the fluid and the solid part), v is the filtration velocity, p° is the pres- 
sure and K is the permeability tensor, depending on the pore structure. In real-world applications 
if is a measurable quantity. In our periodic model it can be computed from the cell problem 

—fiAw- 1 + Vtt 3 = cj in y 

divw 3 = in y (8) 
(w 3 , 7T 3 ) is y — periodic 



1 where e stands for the pore size 



2 



by 

Kij = [ w) = n [ Vw'V^ , (9) 
Jy Jy 

where y c]0, 1[™ is the fluid part of the unit cell (the period) Y =]0, 1[" and S is the boundary of 
its solid part A = Y\y of Y. The viscosity n is equal to 770, in case 7 < 1 and to -q^ if 7 > 1 and 
rjoo > 0. Defining 

u {x,y)=Y^w i {y)^{x) , p'fotf) = X>*(y) |£(x) - 

i=l * i— 1 * 

the above filtration law (^) and the cell problem can be written in the coupled form 

-fj, A y u° + Vyp 1 = f{x) - V x p°(x) in Q x y (10) 
div a u° = in fl X y (11) 
(it^p 1 ) is y — periodic in y (12) 



divz / u°dy = in (13) 



u° = onflxS (14) 



n- (^J u°dy^j =0 on 



(15) 



Two functions u° and p 1 are, in fact, first oscillating terms in the asymptotic expansion for the 
velocity and for the pressure, respectively. 

2. If the viscosity obeys the Carreau's law with 7=1, the homogenized law is the coupled Carreau's 
law 

-div a {ry c (e a (u°))e a (u )} + Vyp 1 = f(x) - V xP °(x) in Q x y 

divyu" = in Qxy (16) 



(vPjp 1 ) is y — periodic in y 
div x (^J u° dy^j = in f2 



>y 

= onfixS 



n- (^J u°dy^j = ondn , 

vc(0 = (m - voo)(i + A|ei 2 )' r/2 - 1 + Vo* ■ (17) 

Here p° is the pressure and 

v = I u°(x,y) dy 



Jy 

is the filtration velocity. Functions u (x,y) and p l {x,y) are oscillating in y and, in analogy with 
the linear case, can be seen as the first oscillating terms in the asymptotic expansion of the flow. 

3. In case of power-law, the homogenized law is the coupled power-law 

-div a {r ?p (ej / (u°))e a (u )} + Vyp 1 = f(x) - V x p°{x) in fi x y 
div y u° = in rixy 

{v° ^p 1 ) is y — periodic in y (18) 



div x (^J u° dy^j = in fi 
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i° = on O x S 



n (^J u°dy^j = on dfl , 



v P (o=m r - 2 t- (19) 

In case 1 two scales are separated and the filtration law is completely macroscopic. All the microscopic 
information are contained in the permeability tensor K , which can be measured in applications. In cases 
2 and 3 the situation is different and two scales are still coupled. Even if, mathematically speaking, those 
results solve the homogenization problem for such fluids, to get a physically relevant result one needs to 
decouple the problems and to find the filtration laws in the usual macroscopic form 

v=U{f-V xP °) (20) 

where 

v(x) = / u°(x,y)dy . 

Jy 

Using the idea from pi] , B and JL2| to decouple the homogenized problems we define the auxiliary problem 



-divy{r) a (e y (w£))ey(w£ )} + V V 7T£ = £ in y (21) 

divywj = in y (22) 

(w^,tt^) is y — periodic (23) 

= on S , a = c,p . (24) 

With (w^,tt^) wc define the permeability function U : R" — > R" by 

W(0= / Mv)dy , (25) 

Jy 

and we obtain formally the filtration law (p0|). From the definition we can obviously conclude that the 
function IA is odd, i.e. that 

w(-0 = -w(0 . 

From (|l^) we also get 

div v = in f2 , v ■ n — on 9f2 . (26) 

To show that ( p"8| ) and (pp[)- (p6[ ) are equivalent, we prove in sections ^ and ^ that W is monotone and 
coercive, consequently, that (|20|)-(|2^) is well posed. 

We also perform the qualitative analysis of the permeability function U. In case of coupled Carreau's 
law, we find its Taylor's expansion. 

2 The Carreau's coupled homogenized problem 

In this section we use the notation r\ for the Carreau's viscosity ?y c , i.e. we omit the index c. We begin 
this section with statement of the main results 



3 Statement of the main result 

To prove that our formal separation of the coupled problem is meaningful we have to prove that our 
macroscopic problem has a unique solution. 
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Theorem 1 Let U : R n — > R" 6e defined by @). TTien </ie macroscopic problem 



divU(f -Vp°) = in ft (27) 
n-U(f- Vp°) = on ,9ft (28) 

/ias a unique (up to a constant) solution p° € VF 1 ' 1 " (ft), /or any f £ L r (ft)™. 

In order to approximate the permeability function IA in vicinity of we compute its Taylor's expansion. 
First two terms can be written in the form 

W(0=^+^(%-r?oo)A(2-r) £ H^^ e e m + 0(\^\ 5 ) , (29) 

m.j,k^— 1 

where if is the classical Darcy's permeability and the coefficients H*™ are defined by (|5^). S uch 
approximate filtration law, given by the polynomial permeability 



1 " 
V(0 = KZ + ~(r ]Q -r }oo )\(2-r) ]T e m 

is still well posed: 

Theorem 2 Let f £ £ 4 (ft)™ . The problem 

divV{f-\7q°) = mil (30) 
n ■ V(/ - V(?°) = on 9ft (31) 

/ias a unique (up to a constant) solution q° € VF 1,4 (ft). 

3.1 Taylor's expansion of the permeability function 

We want to find the Taylor's expansion for U in vicinity of 0. Since the function is odd all the derivatives 
of a pair order D a U(0) , \a\ = 2m , m £ N are equal to 0. Deriving (|2l| ) with respect to £j we find 



Theorem 3 Let K be the Darcy's permeability tensor computed from the local problem 

— ryoAw 7 + V7T 3 = ej in y 

divw j = m y (32) 
(w 3 , it 3 ) is y — periodic 

by 

K i:j = [ w) = 7] [ Vu^Vtu*' . (33) 
Jy Jy 

Then 

V € W(0) = X (34) 



To prove that we proceed as in |J , || or |T^] . We begin by recalling the result from |1J] : 
Lemma 1 There exist a constant c\ > such that : 

\e(v — t/*)!^ r 

L!iy) < [rj(e(u))e(u)-ri(e(v))e{v)]e(u-v) (35) 



1 l + |e( U )|^ ) + Hv)\ 2 L -r r {y) ~ Jy 

for any 



u, v £ W = {<f> £ W 1 ' r (y) n ; dive/) = , <p is y - periodic , <\> = on S} 
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Lemma 2 Let (i/^,7r e ) e W x 6e denned by $H\)-(2£). Th 



for any £ sitc/i i/ia£ 



|e(^)U. w < 2^|y| 1 ^ 1^1 , (36) 



if I < £1 



where C r is the Poincare-Korn's constant in W (i.e. such that \4>\L r (y) _• C"*l e ( ( / , )|L' ~(y) V S W,). 
Proof. Multiplying ( pl[ ) by and integrating over 3^ we obtain 

r |e(tef)lir W 



Now (^5|) gives 

c r Mw 6 )\ Lr{y) <cnyr^(i + ie( We )i 2 L 7 ( ^) ici . 

If |e(u^)|L>-(3;) < 1 then ( |§6| ) obviously holds. On the other hand, supposing that ^(u^Il^) > 1, implies 

cJ|eK)U- (y) < 2C|y| 1 -"| e (^)|^ ) |?| 

and 



|eK)|27 ( 1 w <2^|3;| 1 ^|e| 



which contradicts the assumption that 

\z\< 



2C r \y\ 1 ~- 

Proof of theorem 3. To prove the claim we only need to show that, as \h\ — ► 

n 

\w h -W(h)\ w = o(\h\) , W{h) = Y,h i w i (37) 

i=l 

because 

- Jf /i| = | / K - < CK - W(h)\ w ■ (38) 

In fact we shall prove that 

\U{h)-Kh\ = 0{\h\ 3 ) (39) 

which implies, not only that WU(0) = K, but also that V 2 U(Q) — (which we knew before because U is 
an odd function). Using the standard regularity result (see e.g. Q) for the problem ([?2j) we see that 

n 

\W(h)\ H s {y} <J2\h\\w%s (y) <C\h\ . (40) 

i=i 

and therefore 

\W(h)\ cim <C\h\ . (41) 

Functions W(h) and n(ft) = ^ZILi ^« 71-1 satisfy the system 

-div{»Ke[W(ft)]}e[W(ft)]} + VII(/i) = ft + i4 (42) 
divW(fc) = (43) 
(W(ft), n(ft)) is y - periodic , =0 on S , (44) 



(i 



where 

R h = div{( m - Voo )[(l + \\e[W{h)]\ 2 y' 2 - 1 - l]e[W(h)}} . 
An easy computation leads to 

(i + xey' 2 - 1 - 1 = (r/2 - \)t 2 [\i + Tt 2 y' 2 - 2 dT < ct 2 . 

Jo 

For 1 < r < 2 we get using (|d|) 

Rh < C\4>\ w i,r {y) \ \e[W(h)]\ 3 L3r , (y) < C\h\^\ w ,,r {y) . (45) 
Subtracting (|2lJ) from ([42]), multiplying by — Wh and integrating over ^ we obtain 

J= [ {v{e[W(h)]}e[W(h)]-r]{e(w h )}e(w h )}e(W(h)-w h ) = 
R h (W(h)-w h ) <C\h\ 3 \W(h)-w h \ w . 
An application of lemma 1 gives 

_ r \W(h)-W h \ 2 w 

- 1 i + c\h\ ■* 

As we have seen the second derivative in £ = is 0. Third can be computed from the auxiliary problem 

-jjoAw 1 ^ + Vir £jk = (46) 

= fad - J7oo)(r - 2)Adiv{e(w fe ) • e^) e(u/) + (47) 

+e(w k ) ■ e(w e ) e(w j ) + e(w e ) ■ e(w j ) e(w k )} in y (48) 

divw ijk =0 in y (49) 

ttW*) i s y _ periodic (50) 



by 



Since 



d 3 U(0) 

d 3 U(0) m 



r/o / Vw ijk Vw m = (51) 



= (rjo - Voo )(2 - r)X [ {e(w k ) ■ e{vP) e(w e ) ■ e(w m ) + (52) 

Jy 

+e(w k ) ■ e{w l ) e(w j ) ■ e(w m ) + e(w l ) ■ e(w>) e{w k ) ■ e{w m )} = (53) 

= (vo Voo)(2 r)\(Hfc + + H k p , (54) 



where 



H% = Hjjf = Hi? =Hit= <w k ) ■ e{vP) e(w e ) ■ e(w m ) . (55) 

Jy 

This gives the expansion for U in the form 

U(0=K^+l(f l0 ~ Voo )X(2-r) H *kZAk&e m + 0(\S\ 5 ). (56) 



m,j,fc/-l 
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Proposition 1 Let H?™ be defined by Then holds. 

Proof. We use the same idea as in the proof of theorem We define 

n n 

i+l i.j,k—l 
n 1 n 

Q(/l)=^7T < /l < + - ^ TT^hihjhk . 



2 

i+l i,j,k=l 



Now 



-diT{i)r(e[7(A))]e(V(li))} + VQ(/i) = 

= ft - (m - ??oo)div{[i + AKn/O)! 2 )'' 72 " 1 - l]e[W]} + 

A ™ 
+-(r - 2) (770 - ?7oo)div{ ^ e{w l ) ■ e(w j ) e{w m ) hihjh m } . 

It only remains to estimate 

J = (m - ??oo)div{[l + MeiVih))^)^ 2 - 1 - l]e[V(h)}} + 

A ™ 
+ -(r-2)(r ? o-r ?00 )div{ ^ e(w l ) • e(w J ) e{w rn ) hih 3 h m } . 

A simple Taylor's formula gives 

(i+x\er /2 - l -i = lm 2 +o(xv) . 

After recalling that e{w l ) £ C(y) nxn , e(iy ijfc ) G C(3 ; ) n>< ™, we get for any y £ y 

(m ^)div{[i + A| e (v(/o)| a ) r/2_1 - l]e[y(/i)]} = 

A ™ 
= -(r-2)( ?7o -7 7oo )div{ ^ eK)-eK) e ( W m )^/ 1 ^„ l } + 0(|/i| 4 ) 

Arguing as in the proof of theorem ^, we get the claim. £ 

Remark 1 It can be seen from the definition that Hjj > , Hf^ > 0. Furthermore, for 

n p n n n 



i,j,k,e=l y k=l 



Z(h)-h = / |Ve(»%| 2 >0 . 
k=i 

To compute further order terms we proceed in the same way and we get another auxiliary problem 

-r) Aw ijHp + Vir ijMp = (57) 

= fad - ??oo)A(r - 2)div{ (58) 

e{w l ) ■ e(w j )e(w Hm ) + e{w k ) ■ e{w l )e(w^ m ) + e(w l ) ■ e(w e )e(w jkm ) + (59) 

e{w l ) ■ e(w p )e(w lke ) + e(w j ) ■ e(w k )e(w Um ) + e(w J ) • e{w l )e{w lkm ) + (60) 

e(w j ) ■ e{w p )e(w M ) + e(w k ) ■ e{w l )e{w ijm ) + e(w k ) ■ e(w p )e{w iji ) + (61) 



S 



e{w e ) ■ e(w p )e(w ijk ) + 



(62) 
(63) 
(64) 
(65) 
(66) 
(67) 
(68) 



(e(w jkt ) ■ e{w l ) + e{w lM ) ■ e{w J ) + e(w ije ) ■ e(w k ) + e{w ljk ) ■ e(w e ))e(w p ) + 
(e(w jkm ) ■ e(w l ) + e{w ikm ) ■ e(w J ) + e(w ijm ) ■ e(w k ) + e(w^ k ) ■ e(w p ))e(w e ) + 
(e(w 3lm ) ■ e{w l ) + e{w ltm ) ■ e{w j ) + e(w ijm ) ■ e{w l ) + e{w iji ) ■ e(w p ))e(w k ) + 
(e(w Mm ) ■ e(w l ) + e(w tlm ) ■ e(w k ) + e(w tkm ) ■ e(w e ) + e(w iU ) ■ e{w p ))e{w 1 ) + 
(e(w Mm ) ■ e(w 3 ) + e{w llm ) ■ e(w k ) + e{w lkm ) ■ e(w e ) + e(w jki ) ■ e{w p ))e{w t ) + 
A(r-4){ 



{e{w t ) ■ e(w k )e(w j ) ■ e(w i ) + e(w e ) ■ e(w j )e(w k ) ■ e(w i ) + e(w e ) ■ e(w i )e(w k ) ■ e(w j )}e(w p ) + (69) 
{e(w p ) ■ e(w k )e(w j ) ■ e(u/) + e(w p ) ■ e(w j )e(w k ) ■ e{w i ) + e(w p ) ■ e(w i )e(w k ) ■ e(w j )}e(w e ) +(70) 
{e(w p ) ■ e(w e )e(w j ) ■ e{w i ) + e(w p ) ■ e(w j )e(w e ) ■ e(w i ) + e(w p ) ■ e{w i )e{w e -) • e(w j )}e(w k ) + (71) 
{e(w p ) ■ e(w e )e(w k ) ■ e(w l ) + e(w p ) ■ e(w k )e{w l ) ■ e(w l ) + e(w p ) ■ e(w l )e(w e ) ■ e{w k )}e(w 3 ) + (72) 
{e(w p ) ■ e(w e )e(w k ) ■ e(w j ) + e(w p ) ■ e(w k )e(w e ) ■ e{w 3 ) + e(w p ) ■ e(w 3 )e(w £ ) ■ e(w k )}e(w l )}}(73) 
divw ijMp =0 in y (74) 
(iu« M P,7r yMp ) is y - periodic . (75) 



Now 




and we get an expansion in the form 



W(0 = + -(vo - ?7oo)A(2 -r) ]T H^SkSt *m + 



where 



jke ~ 



mpq 



Jy 




This computation can be proceeded. 



9 



3.2 Proofs of existence theorems 

To prove the theorem ^ we first prove the following lemma that gives coercivity and boundness of the 
differential operator 

A{4>) = -dbrU(f -V<f>) ■ 



Lemma 3 There exist constants M, C, c > such that 

U(0 ■ t > clel"' for any \£\ > M (76) 
|W(OI < Cl^r'" 1 for any ? £ R" . (77) 

Proof. The upper bound is easy to prove, since we have, due to (|35|) 

r |e(^)|ir W 



1 



<c\t\ 



To prove the lower bound ( [r6|) we define — l^ 1 r and — |£| 1 £. Those functions satisfy the 
system 

-div y {n ( {e y (v 6 ))e y (v 6 )} + V y q ( = in y (78) 

divj,^ = in y (79) 
(«{,?{) is y — periodic (80) 
Df = oeS , (81) 

where 

Applying again ( |35| ) we obtain 

|eK)U-(y) < r 

\Z\W-"W-r) + \a(Vt)\ L r (i y ) ~ 

so that we can extract a subsequence such that — ► oo , |Cn| _1 £n ~~ * Coo > |£oo| = 1 and 

V£ n — 1 weakly in W as n — > oo . 

Using the Minty's lemma we get that Voo satisfies the system 

-(no - Voc)y /2 _1 div y {|e y (u oc ))| r ~ 2 e a (w oc )} + Vj,^ = Coo in y 
diVyWoo = in y 
(woo,9oo) is y- periodic 
Woo = on 5 , 

For any £oo, such that \£oo\ = 1, we can define a continuous function 

Zftoo) = / «oo ^ . 

Furthermore 



X(Coo) ■ Coo = / (tto - Voo) X r/2 - 1 |e(«oo)| r 

and there exists 



y 



Coo = inf T(t) ■ t > 

|r|=l 
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Two functionals 

Jy 

*<»(») = I (vo-Voo)^ 2 - 1 |e(«)|* 



are obviously convex and, due to the Lebesgues dominated convergence theorem, 

- $ao as |£| ->■ oo . 
Since $oo is weak lower semicontinuous in W , we obtain 

lim inf $ e (u e ) = lim inf-^^) + - $oo(fc)]} > $«>(««>) > > . 

|f|-«x) l£|->°° 

But 

$ c K) = ier r '^(e)-c 

which proves the claim. X 

Proof of theorem 3. Lemma [T^ implies coercivity and boundness of the operator A : W 1 ' 7 ' (0)/R 
(W rl,r (Q))' . Its strict monotonicity follows from inequality 



[U{£) - U{t)] ■ {i - t) = I {ri c (e(w 6 ))e(w^ - Vc (w T ))} ■ (w s - w T ) 

\e(w£ - w T )\ 2 Lr(y) 



\<w{]\ 2 Lr r {y) + |eK)lirfy) + 1 



Now the classical Browder-Minty's theorem proves the existence of the solution, i.e. the surjecivity of 
our operator. As it is strictly monotone, it is also injective, i.e. the solution is unique. X 

If we approximate IA by taking only first two terms in its Taylor's series, i.e., if we consider 

1 - 



instead of U, we get the problem fl30|)-(|3l|). 

For such approximation of the filtration law we have theorem ^, from the beginning of this section. We 
are now able to prove it. 

Proof of theorem ||]. We define the formal differential operator 

BfeO=divV(/-Vp) . 
Since V is a polynomial of order 3, it satisfies the estimate 

|V(£)|<C(|£| + |£| 3 ) . 
Consequently B : W 1A (Q)/R~> (W 1,4 (f2)/R)' is a bounded operator. Moreover 

\BQp)\( W i,*(si)/B.y <C\f - Vv?|i 4(n) . 
Next we prove the coercivity. Let 

n 

E(0 =$>(«/*) & . 
i=l 

Then, for any £ € R™, 

V(0 ■£=/(% \E(0\ 2 + (m - Voo) A(l - r/2)\E{0\ A ) ■ 
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The functions 

5(0= [ \m)\ 2 , no= f m)\ 4 

Jy Jy 
are continuous and strictly positive for £ =^ 0. Therefore there exist 

/si = inf x 5(0 , k 2 = iMWCO ■ 

But than 

V(0 • £ > «i |£| 2 + (770 - r?oo) A (1 — r/2) k |£| 4 , 

implying the coercivity of B on W /1 ' 4 (f7)/R. It remains to prove the monotonicity of B. Then the result 
will follow from the classical Browder-Minty's theorem. 
For any £, r € R" the derivative -DV(£) satisfies 

DV(0 T.T = Jfr.r+ l -\{2 - r)(t» - ?7oo) / [(£(£) ■ £(r)) 2 + 



>' 



-£(£) 2 £(r) 2 ] >]fr-r> fc |r| ; 



The above inequality proves that the equation ( p0| ) is uniformly elliptic , i.e. that the operator B is 
monotone. X 

The approximation V of the permeability function is valid either for small |£| or for small A. For large A 
and the Carreau's auxiliary problem (pT|)-(p4"|), with r] c , is close to the power-law auxiliary problem 
(pl|)-(|24|), with r/ p and /i = (r/o — T] 00 )X r '^ : ^ . We study the power-law auxiliary problem in the next 
section. 

We refer to section |6.3| for numerical computation of the permeability function U and its approximation 
V and their comparison for sufficiently small 



4 The power-law coupled homogenized problem 

In this section we also start with statement of its main results and we adopt the notation r\ for the 
power-law viscosity r\ v . 

5 Statement of the main result 

From j2l|)-([25|) we get by substitution that the permeability function is homogeneous in the sense that 

H(A£) = |A| r '- 2 AW(0 , AeR . (82) 

implying 

HCr'- 1 < \U(i)\ < Mltf'- 1 (83) 

with 

m= inf \U(0\ , M= sup \U{0\ ■ 
I€l=i |£|=i 

Now, as in the previous chapter, ( |S3"| ) and (§3) imply that: 
Theorem 4 For any / £ L r (O)™ i/ie macroscopic problem 

divU{f - Vp°) = m (84) 
n-W(/ - Vp°) = onffi (85) 

/ias a unique (up to a constant) solution p° € W 1,r (f2). 
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The standard engineering model (see e.g. 0, |l5[]) suggests that the permeability function U can be 
written in the form 

U(0i = \{KZ) i \ r '- 2 (K0i , i = l,...,n , (86) 
where K is a permeability tensor. 

To check that conjecture we define the conjugated function G by 
such that 

Now, if the conjecture ( ^6| ) is true, the function G is linear. Homogeneity ( j82| ) of U implies that G is 
homogeneous with order 1, i.e. that 

G{\0 = AG(0 , AeR . 

Function homogeneous with order 1 is linear if and only if it is differentiate for £ = 0. 
Function G is differentiable for any £ ^ 0, but, in general, not for £ = 0. That means that the conjecture 
( |86| ) is, in general, false. In case of unidirectional flow it is differentiable and, therefore, linear. In 
case of a porous medium consisting of a net of narrow channels, G allows an asymptotic expansion in 
powers of the thickness of the channels, with the first term being a linear function. Numerical examples 
show that G is linear (or differentiable for £ = 0) in case of periodic porous medium with the solid part 
being an array of rectangular obstacles (i.e. with the fluid part being a net of perpendicular channels). 
On the other hand, in case of spherical or elliptical obstacle, the function G is not linear and conjecture 
( |86| ) does not hold. 

5.1 Study of the permeability function 

We next recall the result from juj, that will be used in the sequel: 
Lemma 4 There exist a constant c\ > such that : 

\e(v-u)\ Lr{y) < f l n ^ u )) e (u)-r){e{v))e{v)]e{u-v) (87) 



^eiut^ + Hv)] 2 -^ ~Jy 
for any 

u, v <G W — {(j) e W 1 - r (y) n ; divcf) = , 4> is y - periodic , (j) = on S} 
It implies that, for £, r 6 R™, £ ^ r, 

\e(w(:-w T )\ 2 Lr ( y) 



[W(0 - U{r)\ ■ (£ - r) > c ^J^ %_ r > 



Before we continue we give two simple examples mentioned in the introduction of this section. 

Example 1 As it was noticed in ^/, J7j|/ , /|]/, if the flow is unidirectional, e.g. in the direction e\, then 
W(0 =J(6) ei, with 

J(di) = \Zi\ r '- 2 Zi • 



This result corresponds to the usual models that can be found in the engineering literature (see e.g. 15]. 
• 
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Example 2 In case of a porous medium consisting of a network of thin pipes we can find the asymptotic 
form of our function li and, again, it corresponds to the usual engineering model from ffij and Wjj. 
Suppose that n = 2 and that y = y$ consists of two narrow perpendicular channels 

Pi =]0,l[x ] -6/2,6/2 [ , Pi=] -6/2,6/2[x]0,l{ , 

with thickness S. We consider our auxiliary problem in such unit cell ys = P} U P$ 

-dw y {r) p (e y (w^))e y (w s ^} + Vj,tt| = f in y s (89) 

div y wl = in y s (90) 

(u>£,7Tf) is y$ — periodic (91) 

w 6 = onS 5 , (92) 



where S = d(Pg U P$). We study the limit as 5^0. The similar problem was considered in f/p\ j for 
the Newtonian fluid and an explicit formula for the permeability was found. We are going to do the same 
thing with our nonlinear permeability function ti s {^) — JyW^. As in /pi]/, we seek an approximation of 
(wi,iri) in the form 



% 2 /<5) lar'^i ei m Pi f y 2 & in Pf 

9(yi/6) |6| r '- 2 C 2 e 2 »P 2 S ' 1 ' \ yiti m P & 2 



where 

In each P\ we have the Poincare-Korn's inequality 

\w\ - w{6)\ L r {P i ) < C8\e{w\ - w{S))\ L r (Pl) 

and the trace inequality 

H - wWlLrpt) < C 5 1/r ' |e(tfl| - w{8))\ L r {Pi) , 

where 

By a direct computation we obtain 

A* [__{\e(wl)\ r - 2 e(wD - \e( W (S))r 2 e(w(S))}e(wl - w(8)) = 
{n(5)n - ti\e(w{8))\ r - 2 e{w(8))n} (iu| - w(5)) . 

We use ([8lj ) to estimate from below the left-hand side and the above inequalities to estimate from above 
the right-hand side. It leads to the estimate 

U s (0 - (9) J(0 \<C8 1 ' r , 

where 

{0) = f 2 6(s)ds = —l—(V2/^y'-i , 



and 

^(C) = (l6l r '- 2 a,l6l r '- 2 6) 
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Therefore 

mOi = s r,+i —^—(V2/ny'- 2 le.r'- 2 & + o{s r ' +i+i ' r ) , * = 1,2 

(r + 1) fj, 

i.e. 

Wa(0 « . 



We are now tempted to believe that the permeability function can be written in the form (p6|). 

First of all, we can not hope to have K that depends only on the geometry of the porous medium, as in 

the linear case, because in two above examples it depends on the flow index r. 

We know that this conjecture is true in case of unidirectional flow and, up to some order of approximation, 
for a network of narrow channels. 



To verify whether ( p6[ ) is true or not, we define the conjugated permeability function 

G{0i = \m)i\ r - 2 U{0i > i = l,...,n . 

As r — 2 < and U(0) = , G(0) is not correctly defined by the above formula. However, since U satisfies 
@, and @we have 

m\i\ < \G(0\ < M\C\ . 

Therefore we define G in £ = by G(0) = and we get the function G G C(R) n . Function G is obviously 
homogeneous 



we refer to the section 5.4 where a numerical example shows that in case of bundle of perpendicular, 
not necessarily narrow, channels G is differentiable and linear, i.e. that (|8^) can be expected to hold 
in rectangular geometry. That can be explained by the fact that we have n unidirectional flows with, 
almost, no interaction between them. In fact, example ^| proves that the interaction between two flows 
exists, but it has some lower order. 

As it was noticed in |l3| the function U is differentiable for any £ ^ 0. We are, therefore, in situation 
very close to (|S6|). Moreover 



wl(y)}idy , (93) 
y 



where 



+ \e v (w s )r 2 e v (wl)} + Vyirl = e, in y (95) 



(we , 7r|) is y — periodic 

wl = on S . (96) 



- M div,{(r - 2)\e y ( Wi )\ r -* [e(w t ) ■ e(«#] e y (w 6 ) + (94) 
divyi/;^ =0 in y 

4,4) 

Since this problem is linear it is easy to prove the following result: 

Proposition 2 For £ 7^ the problem fi9$l)-ft9(\) has a unique solution «;| G V^, where 

V C = {<t> G W^ r (y) n ; leK)!*-^) G L 2 (tt) , ^ = } 

and 

W^ r (n) = {(/) G W-^CV) ; ^ is y ~ periodic , w^OonS}. 
Furthermore satisfies the a priori estimate 

\e(wi)\ Lr{y) <C\e(w c )\ 2 - r (y) . (97) 
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Proof. We first notice that the definition of has a sense since, due to the e(iyg) ^ (a.e) in y. 
Now is the Banach space equipped by the norm 

l<% = |efa)lz-(y) + I \e(w 6 )\i- 1 e(4>)\ L ^ (y) . 
The operator T :Vjt VI, defined by 

T4> = -fidiv y {{r - 2)|e y (w c )| r ~ 4 [e(w^ • e(0)] e y (w ( ) + \e y {w(:)\ r - 2 ey{<i))} 
is coercive since, using the Holder's inequality, we get 



(TM >(r-l) / |eK)r^|e(0)| 
r- 1 



> 



{| |e(«7 C )|*- 1 e(0)|i» c „ 3 + (Hw^y^-^)^} 



2 

Now T € £(Vj, W ) is linear and coercive. Therefore it is bijective. The estimate ( |97| ) follows from 

Lemma 5 TTie matrix Vj£Y(£) , £ ^ is symmetric. 

Proof. Multiplying (^3) by W and integrating over 3^ we obtain 



^(0 = / {(r- 2)|e y (^)p 4 [e(w ( ) ■ e(w\)] [e y (w^ ■ e y {w{)} + 



dU 

+\e y (ws)\ r -\(wl)-e y (wi)} = ^(0 .* 



>' 



Now U is differentiable, its derivative is defined by (|95[)-(|96|) in any point £ g R™ except for £ = and it 
is symmetric. Consequently G is differentiable for any £ 7^ with 

^(0 = (r-l)|M(O,f- a ||(0 . (98) 

Due to the estimate 

\U{j)-U{Q)\ ,-2 

we see that V^W(O) = 0. It is, therefore, not possible to linearize U in vicinity of £ = 0, as we did 
in the Carreau's case. Once again we have the problem with G for £ = 0, since ([)8j) becomes an 
undetermined expression of the form 2. But, we obviously have 

|V«G(0| <C , 

so that G € C'^R"). 

Unfortunately, this is as far as we go with (|8^ ) . The two numerical examples found in section |6.4| show 
that, in general, G is nonlinear. That, of course, means that it is not differentiable for £ = 0, which can 
also be seen from those numerical computations. 



6 Numerical Experiments 



To illustrate the various theoretical results of this article, we performe in this section numerical compu- 
tations. After a short description of the numerical schemes used to solve the linear and nonlinear Stokes 



problem, we describe the different pores geometries considered. The following section (6.3) is devoted to 
the numerical experiment in case of Careau's viscosity, and finally the last section (6.4) is concerning the 
power-law viscosity. 
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6.1 Numerical methods used 



We have to compute the solution of the nonlinear Stokes system (pi|-p4|) and for the case of the Carreau law 
viscosity, the solution of linear Stokes problems. Two difficulties have to be solved: the incomprcssibility 
condition and the nonlinearity in case of nonlinear viscosity law. An Augmented Lagrangian algorithm is 
used to satisfy the incompressibility condition div u — 0. We use a finite element method to discretize the 
Stokes problems: a P 2 ° is used for the velocity field, and a P\ discontinuous element is used for the pressure 
field. Those elements are known to give good results with Augmented Lagrangian algorithm. Finnaly, 
the nonlinearities caused by the viscosity law are solved by a Newton Raphson procedure. We refer to || 
for the description of the global algorithm. The porous cells are discretized by a mesh containing 1000 
triangles. 

6.2 The Porous media considered 

We consider three models of porous media to make numerical experiments: Each of them arc characterized 
by the solid inclusion in the periodic unitary cell. As shown in figure ([!]) , we consider three geometries of 
the inclusion: a cylindrical inclusion, a square inclusion, and finally, an ellipsoidal inclusion which gives 
an anisotropic porous medium. In the following, we refer to GEOM1 for the prous medium with circular 
inclusion, GEOM2 for this one with square inclusion and finally, GEOM3 for the last porous medium 
with ellipsoidal inclusion. 




Figure 1: Periodic unitary cell: GEOM1 on the left, GEOM2 in the middle and finally GEOM3 on 
the right of the figure. 



6.3 Numerical Simulations for Carreau's viscosity law 

In this section, we propose to solve (^2l|-|24|) for a sampling of external force £ and its Taylor's expansion 
around the origin when viscosity is governed by a Carreau's law. At first, we describe parameters of 
the Carreau law and the numerical experiments procedure. In applications (see e.g. |]l|) 7700 is either 
very small compared to 770 (polymer solutions) or even equal to zero (polymer melts). For this reason, 
neglect r]oo and consider a simplified Carreau law viscosity with rjoc = 0. We chose arbitrarely 770 = 1 
and consider two values of A: Ai = 1 and A 2 = 100 to take into account the different comportment for a 
viscosity law: one near a constant law, and the other one far of the constant case. Finally, we consider a 
flow index r — 1.5: 

rj(e(u)) = (1 + A.le^)! 2 )"- 25 with Ai = 1. or A 2 = 100. 

Because theoretical results are given for small £, we consider at first, the solution of the nonlinear problem 
- |24| ) for |£| < 1. Numerical solution are performed for a sampling of £ defined by: 

f# = *(l C os— , -sin—), l<i<n = 5, < j < m = 8 
n 4m n 4m 

extended by symmetry for m < j < 4m . For each value of £ lJ we compute the gap (resp.Af j ; Af^) 
between the solution Uij of (f2l"|-p4]) and the Taylor's expansion of order 1 (resp 3, 5) where A*- are defined 
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On figure (||) we plot different A^- for fc = 1, 2, 3 for A 2 = 100 and for the first two geometries GEOM1 
and GEOM2. Concerning the first geometry, one may see that we need the third term in the Taylor's 
expansion to obtain good results. But, in the case of big square inclusion (GEOM2), for this range of 
|£|, the global filtration law is very close to the Darcy's law: there is no real gap between Darcy's law 
and the Taylor's expansion of any of three orders. To underline those comportments, we have drawn on 




Figure 2: Gap between U and its Taylor's expansions of order 1, 3 and 5 for GEOM1 on the left, and 
for GEOM2 on the right (for A 2 = 100). 



figures (f|{|) the same quantities (eq. (J9g| — 101 ) ) for a fixed direction of £ and for length < |£| < 1. Let 
consider first figure (||) concerning GEOM1. In the case of Ai = 1 one may see that for the range of 
considered £, the filtration law is very close to Darcy's law: the relatives errors, independently of order, are 
smaller than 1%. But in the case of A 2 = 100, one needs the fifth approximation to obtain error smaller 
than 1%. Secondly, when considering figure ([|) concerning GEOM2, one may see that independently 
of the value of A, the behaviour of IA remains almost linear: for Ai = 1 error is smaller than 0.8% and 
for Ai = 100 error is smaller than 8%. To close this numerical study of theoretical result s co ncerning 
Carreau's viscosity, we have computed solution of (|2l|-|2^) and the different gaps (eq. ([99|- |l0l| )) for |£| 
larger than 1, to obtain numerical limit of validity of the Taylor's expansion, for the three geometries 
but only for one fixed direction of £: Z(ei,£) = 60°. Those results are plotted on figure ||. One may see 
that for the three geometries, Taylor's expansions give good results for Ai = 1 for the considered range 
of £, but in the case of larger parameter A 2 = 100 one loses precision of the Taylor's expansion. More 
precisely, the best results are for the third order Taylor's expansion, but the fifth order expansion gives 
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0.20.30.40.50.60.70.80.9 1 0.20.30.40.50.60.70.80.9 1 



Figure 3: Gap between U and its Taylor's expansion for GEOMl:|£| < 1, (ei,£) = 45°, A = 1 on the 
left, A = 100 on the right. 




0.20.30.40.50.60.70.80.9 1 0.20.30.40.50.60.70.80.9 1 



Figure 4: Gap between U and its Taylor's expansion for GEOM2:|£| < 1, (ei,£) = 45°, A = 1 on the 
left, A = 100 on the right. 

very bad results ( the rest of the Taylor expansion of order 7 becomes non negligible) . This suggests that 
for large |£| and A the Taylor's series diverge. 

6.4 Numerical simulation for Power viscosity law 

In this section, we propose to solve (^l|-^4|) for a sampling of £ when viscosity is governed by a Power 
law to illustrate the theoretical results of section [|. At first, we describe parameters of the Power law 
and the numerical experiments procedure. Due to the homogeneity property of the power law viscosity, 
we just consider a viscositiy law with A = 1. We chose arbitrarely a flow index r = 1.5 : 

rj(e(u)) = \e(u)\- 5 

For the same reason of homogeneity, we consider only the solution of the nonlinear problem (pl|-p4[) for 
|£| = 1. Numerical solution are performed for a sampling £■? of £ defined by: 

^ =t(cos £' sin £ ) '°- j - m=8 (102) 

and any extended by symmetry as in the Carreau's case. In two following sections, we will numerically 
show the different behaviour of the filtration law depending on the geometry of porous medium. 
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Figure 5: Gap between U and its Taylor's development of order 1, 3 and 5 for |£| < 10, for GEOM1 on 
the left, and for GEOM2 in the middle and GEOM3 on the right. 

6.4.1 Numerical results 

This section is devoted to numerical results in the case of power law to exhibit the linear or nonlinear 
comportment of the function G(£) defined by its components G(£)j = \Ui{£)\ r ~ 2 Ui{£) . On the figure (0) 
we plot the first component Gi of the function G for the three geometries £ (due to the symmetry of 
inclusion, the comportment of the second component of G is deduced from the first one by an appropriate 
symmetry) . It may be noticed that the comportment of Gi as a function of seems linear for the case 
of the square inclusion (GEOM2), but one sees that in the case of the spherical or ellipsoidal inclusion 
( GEOM1 and GEOM3) , this is not the case. To underline this effect, we have computed the slope A 




Figure 6: First component Gi(£) of the function G for GEOM1 on left GEOM2 in the middle and 
GEOM3 on right 

of the straight line y = A£i which approaches, in the least square sense, the function Gi for the each of 
the three geometries. This slope is solution of the following problem: 

N 
7 = 1 

and is given by: 

N 



7 = 1 
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Then we have computed the relative distance A between the straight line and the graph of Gi : 

In table)!] we give the values of A and A: One sees that A is the smallest for the square inclusion 





GEOMi 


GEOhh 


GEOM 3 


A 


-0.24083 lO" 1 


-0.63016 10~ 2 


-0.41984 10" 1 


A 


0.99781 lO" 1 


0.44640 10" 1 


0.27976 



Table 1: Slope A of the straight line and relative distance A between this line and G\ for the three 
geometries 

(GEOM2), (of order 4%), and is no negligible for the other geometries (10% for the spherical inclusion 
and 27% for the ellipsoidal inclusion). Finally, we have drawn on figure (0), the pointwise distance 
between the straight line y = A£i and the function G(£)i = |Wi(£)| r ~ 2 Wi(£) for the three geometries. 
One may see that the gap between the straight line and the function G\ is maximal at the point where 
U\ is changing sign (for angle 7r/2 for GEOMI and GEOM2). This effect may be understood as the 
non-differentiability of hi in the origin. 
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Figure 7: Pointwise gap between G\ and its least square approximation y = A£i and Gi(£) for 
GEOMI, GEOM2 andGEOM3. 
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